38 38 X& 第 4 期 E h GE dE Vol.38, No. 4 
1995 年 11 H ACTA ENTOMOLOGICA SINICA Nov., 1995 


重 数 量 与 宿主 数量 和 气象 因子 的 天 系 
地 k à 


《北京 师范 大 学 数学 系 北京 100875) 


张 万 R 


《内 蒙古 伊 克昭 盟 地 方 病 防治 站 KEE 017000) 


马 Y 名 


《吉林 省 地 方 病 第 一 防治 研究 所 ”白城 1370005 


摘要 ”根据 内 蒙古 自治 区 鄂 托 克 旗 和 鄂 托 克 前 旗 1975 一 1989 年 长 爪 沙 鼠 密度 、 重 指数 监测 数 
据 和 本 地 区 气象 站 的 7 项 气象 因子 资料 ， 分 别 求 出 了 蚤 指数 与 鼠 密 度 的 直线 和 曲线 的 回归 模 
型 ,与 气象 因子 的 最 优 回 归 子 集 模型 和 标准 回归 模型 ,给 出 了 鼠 孟 因子 和 气象 因 于 图 的 典型 相 
XA. He. 宿主 数量 变化 导致 重 指数 变化 ;气象 因子 综合 影响 蚤 指数 ;相对 湿度 和 地 诡 温 
度 是 影响 蚤 数量 变动 的 重要 因子 ;气象 因子 对 重 指 数 的 影响 大 于 对 鼠 密 度 的 影响 。 


关键 词 际 指数 ,长 爪 沙 巾 ,气象 因子 ,数学 模型 


重 数 量 与 宿主 数量 ， 已 做 过 一 些 研究 。Cpnpanos P. T. (1963) YIRE TAKEM 
Rhombomys opimus 数量 与 其 寄生 重 数 量 的 调查 中。 马 立 名 《〈1988) 研究 过 达 乌 尔 黄 限 
Citellus dauricus BRAZ Cricetulus barabensis, Xi Hj UTE Marmora himalayana 
WU RH Ochotona daurica 数量 与 蚤 指数 关系 只。 李 书 宝 (1988) 给 出 了 达 乌 尔 黄 
鼠 密度 与 蚤 指数 的 相关 分 析 名 。 李 体 来 等 (19?3) 研 究 了 松 辽 平原 达 乌 尔 黄 鼠 密 度 和 气象 
因子 与 蚤 指数 关系 中。 本 文 将 从 多 种 角度 ,对 长 爪 沙 鼠 Meriones unquiculatus BE, 3 
旨 数 和 气象 因子 的 关系 进行 研究 。 

1 材料 与 方法 

样 地 选 在 内 蒙古 自治 区 鄂 托 克 旗 查 布 苏 木 和 哪 托 克 前 旗 布 拉 格 苏 木 , 地 处 东经 
106°42 一 107°45', 北 纬 38°10' 一 39。 根 据 土壤 、 地 形 ， 植 被 将 调查 区 划分 为 4 种 栖息 地 
类 型 。(1) IEE Stipa breviflora, TĦ Cleisteogenes songorica 为 代表 的 沙砾 
质 高 平原 台地 。(2) 油 萝 4rteraisie ordosica IRRF RW G) AIR Kalidium 


gracile, HH] Nitraria sibirica, EE Ei Achnatherum splendens WERRIET Ho 
(4) 藏 锦 鸡 儿 Caragana tibetica 合 地 。 


本 文 于 1993 年 8 Hk ie 
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在 4 种 栖息 地 内 ，1975--1989 年 每 年 春季 4 一 5 H, KE 10 一 11 月 进行 长 爪 沙 鼠 数 
量 调查 ( 缺 1976 年 .1978 年 调查 资料 )。 样 方 以 公顷 为 单位 ,采用 24h 与 形 夹 法 ， 对 捕获 
只 数 进行 登记 。 每 年 两 次 调查 数据 的 均值 作为 鼠 密 度 y 寄生 重 数 量 用 指数 y, 表示 。 气 
象 数 据 取 上 自 鄂 托 克 前 旗 吉 拉 气 象 站 。 气 象 因 子 取 ; x, = 年 均 气 温 、z 一 年 均 相 对 湿度 、 
v, — FARKE, n 一 年 均 气 压 ,x; — 年 均 地 表 温 度 (0cm)、x6 一 年 燕 发 量 、z 一 年 日 
照 。 监 测 区 在 气象 站 半径 为 50km 的 范围 内 。 数 据 见 表 lo 
表 1 1975—1989 年 长 爪 沙 鼠 鼠 密 度 和 气象 资料 


年 份 | Res 得 指数 ”平均 气温 ”平均 相对 ME 水 平均 气压 平均 地 面 KRE HA 
(只 /ha) 《 匹 ) CO) REC) (mm) (bar) BECC) (mm) Ch) 





1975 59.75 1.98 7.6 56 195.7 864.0 10.1 2517.6 2919.8 
1977 8.08 0.48 7.3 53 294.0 — 864.3 9.9 2362.0 2958.4 
1979 12.76 0.35 7.6 49 267.3 864.0 10.2 2537.8 2868.7 
1980 12.77 0.37 7.4 47 147.3 864.2 10.4 2769.4 3089.0 
1981 0.84 0.42 7.2 50 218.1 864.7 10.1 2353.5 2654.2 
1982 0.30 1.86 7.8 49 147.4 . 864.4 11.0 2611.5 2679.2 
1983 2.40 0.73 7.4 51 198.2 864.7 10.0 2434.8 2472.3 
1984 21.40 0.28 6.4 5] 355.9 863.8 9.3 2179.6 . 2802.0 
1985 118.00 0.27 7.1 51 417.2 . 863.8 10.1 2240.9 2855.5 
1986 30.52 0.74 6.6 50 203.0 867.8 9.7 2388.3 3067.6 
1987 57.30 1.68 8.4 45 181.5 867.3 11.1 2910.5 3103.7 
1988 76.20 1.17 7.0 50 287.4 867.9 10.1 2357.4 2917.2 
1989 162.00 4.80 7.3 57 350.4 867.6 10.7 2350.1 2874.0 
x 43.25 1.16 7.3 51 251.0 — 865.3 10.2 2462.6 2866.3 
s 50.24 1.26 0.5 3 85.4 1.7 0.5 205.5 182.6 
最 小 值 0.30 0.27 6.4 45 147.3 863.8 9.3 2179.6 2472.3 
最 大 值 162.00 4.80 8.4 57 417.2 867.9 11.1 2910.5 3103.7 
cy) . 116.15 107.84 6.9 6 34.0 0.2 4.9 8.4 6.4 


1.1 基本 统计 分 析 
由 表 1 知 鼠 密度 、 重 指数 的 变异 系数 (c. ”) 最 大 , 且 最 高 年 份 鼠 密度 、 蚤 指数 分 别 是 
最 低 年 份 的 540、18 倍 , 表 明 密 度 和 蚤 指数 在 不 同年 份 的 波动 变化 剧烈 。 气 象 因 子 的 e. ys 
较 小 ,年 降水 略 大 。 
12 重 指 数 与 鼠 密 度 的 回归 模型 
首先 , 求 其 直线 回归 方程 为 : 
y; = 0.44552 十 0.01661y, (0) 
显著 性 检验 : F = 8.71 > Fo.0s(1,11) = 4.84 
表 考 虑 曲线 回归 。 在 20 种 曲线 模型 中 马 ， 
选取 下 值 大 的 相应 的 回归 方程 为 
y; = 0.65304 + 0.00012y? (0) 
显著 人 性 检验 : F = 13.81 > F4 (1,11) = 9.65, 
大 我 们 求 二 次 曲线 方程 y; b. 十 by, 十 和 2 好， 并 做 逐步 回归 , 引入 和 剔除 变量 的 临 
3H 7932, AREA). 由 此 表明 ,(2) 恰 为 多 项 式 回归 中 的 逐步 回归 方程 。 
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1.3 ” 重 指 数 与 气象 因子 的 最 优 回 归 子 集 

REER y, 做 为 办 变量 ,气象 因子 ar 为 自 变 量 ， 利 用 回归 分 析 或 逐步 回归 分 析 
方法 ,建立 y 与 x+ 一 *; 的 关系 模型 。 

由 于 7 个 气象 因子 的 一 切 可 能 回归 方程 有 CU 一 1,2,. … ,7) 个 含 ;个 气象 因子 的 
回归 ,共有 27— 12127 个 可 能 的 回归 ,从 中 找 出 一 个 最 好 的 , 它 所 包含 因子 的 回归 方程 即 
为 所 求 , 称 为 最 优 回归 子 集 模型 。 

为 节约 篇 幅 , 仅 写 出 含 i(i 一 1,2, 7) 个 因子 的 最 优 回归 模型 ,结果 见 表 2。 表 2 
是 按 回 归 模 型 中 因子 个 数 的 多 少 依次 排列 的 。 例 如 表 2 中 的 (3), 其 回归 模型 为 

y; = —262.4810 十 0.2762x, 十 0.2695x, + 1.6116z;, 
上 式 是 从 含有 3 个 因子 的 35 个 回归 模 列 中 按 残 差 平方 和 0, 选 出 最 小 的 一 个 所 得 。 其 全 
类 推 。 显 著 性 检验 : (1039 P< 0.10,(2)—(5)28 P < 0.001, (6)—(7)26 P < 0.01, 

JER 2 中 选 出 最 优 的 回归 模型 ,我 们 所 Aitkin(1974) 提出 的 R? 充分 集 的 范围 选 

BUS. i 
Rè = 1] — (1 — RJXYX3!-FpFpn-—p—1i)/n—p-—l1) 


称 满 足 Rz > Ri 的 集 为 充分 集 。 在 R'! 充分 集中 的 因子 才 有 资格 被 选 为 最 优 子 集 。 TAG 
在 满足 该 子 集中 选 R (全 相关 系数 ) 较 大 的 , 且 增 加 溉 子 后 ，R? 增加 幅度 很 小 (或 9 减少 


幅度 很 小 ), 在 此 原则 下 , 选 出 的 即 为 坡 优 回归 模型 
在 本 模型 中 , n—13, p 一 7，R3 一 0.9615 eR Jo ee dod: E: 

1 — (1 一 0.9615)(1 + 7 X 4.88/5) = 0.6985, di 2 中 (2) 一 (7) 可 选 为 最 优 于 集 。 藉 
Ri 较 大 ， "d Ri—R? 增加 幅度 很 小 ( 相 邻 的 RI 之 差 小 于 0.02), 故 (3) 为 最 优 回 归 模 型 。 氢 

AHR H, 

R2 & i(i = 1207.71) 个 气象 因子 的 最 优 回 归 模 型 

， 全 相关 系 REG 

7 AX R* 


(12 b, b, b, b, b, b, b, SKIO 
Ci) — 12.7993 1.36798 0.2940 3.3442 
(2) —32.0276 0.2809 1.8569 0.7855 4.0547 
(3) 262.4810 0.2762 0.2695 1.6116 0.9068 1.7615 
(4) —229.1423  —0.5515 0.2701 0.2301 2.1146 0.9177 1.5552 
(5) —215.8982  —0.9622 20.3031 0.2122 2.1157 0.0014 0.9282 1.3575 
(6) —232.7680 一 1.0936 0.3110 0.2333 2.0650 0.0021 -— 0.0006 0.9331 1.2653 
(7) 一 271.0363  —1.7346 0.3508 0.0059 0.2750 1.7471 0.0066 4.0.0023 0.9615 0.7273 


ZPEGWRIBBUd. 5p» AGAUERCAURIACT UU ERBU2—4.5 中 任 一 盘 , 08181375 4€ 
RAG) EGRI, AXR Rae A [9 C BD 33 27 Da I T 29 352 UG 8L HT, 

将 表 1 EXITRSBUNEUOA EE .SGRUUEUEBUBEUS 

y; = 0.7154x; 十 0.3610x4 十 0.6388 x5 (8) 

其 中 ox; 一 2,4,5) 与 y 分 别 表示 标准 化 内 子 。 
14 和 鼠 重 因子 与 气象 因子 的 典型 相关 分 析 

经 计算 ,得 第 一 组 标准 系数 的 典型 变量 为 ; 

4, = 0.3171y, + 0.7608y, 
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vi = —0.6925x, + 0.8667x, 十 0.6150x, 十 0.3822x, + 0.6842x, 
+ 1.1003x, — 0.2706», 
典型 相关 系数 2, 一 0.9890, 显 著 性 检验 : Qi = 39.75 > XGas (14) 一 36.12, 
第 二 组 标准 系数 的 典型 变量 为 
u, = 1.3006y, — 1.1015y, 
v, = 0.1879x, — 0.3873x, + 0.9585x, 一 0.0057x, — 0.1791 x; 
— 0.1559x, -+ 0.3576x, 
à, 0.8268, T 2E MESA. Q, 8.05 > XL4(6) = 7.840 


2 结果 与 讨论 
21 在 鄂尔多斯 上 暧 温 型 薪 漠 草原, 秀 病 童 蒙 党 亚 种 Nosopsyllus laeviceps kuzenkovi, ln] 
Wæ Xenopsylla conformis, W% Rhadinopsylla dives 分 别 占 主要 宿主 长 爪 沙 
EL: 3E ZR pA 52.3% ,27.2% 10.896 , RIECK, BEC SE ZZ TAURO IR. Aic EOM 
ge GE uk BIKHJo BEEXLARBEESJIAIE. RAES EZ EEMMERI MS PB 
方程 (0) 和 (0) 看 出 ,宿主 数量 的 升降 ， 会 显著 导致 重 指 数 的 升 阵 P < 0.05); BELT EUN 
其 和 伪 指 数 的 单位 不 一 致 , 考 虑 变异 系数 分 别 为 116.15% 和 107.84% ,表明 其 变 守 程度 基 


Zx— 8X , ELSE RES CT EE Ao 
x3 [liées ss 


年 份 3k fa A CUS) [5] J3 48 Ro X 
1975 1.98 1.9113 0.0687 
1977 0.48 0.5738 一 0.0938 
1979 0.35 0.3047 0.0453 
1980 0.37 0.6816 — 0.3116 
1981 0.42 0.3380 0.0820 
1982 1.86 1.6743 0.1857 
1983 0.73 1.0020 —0.2720 
1984 0.28 — 0.2424 0.5224 
1985 0.27 0.5992 — 0.3202 
1986 0.74 0.736] 0.0039 
1987 1.68 1.4274 0.2526 
1988 1.17 1.4048 一 0.2348 
1989 4.80 4.7289 0.0711 


2.2 计算 重 指 数 y, 与 气象 因子 x 一 x; 的 相关 系数 依次 为 rs 一 0.2819, n, s, 一 0.5288, 
ryn, = 0.0536, r,,,, = 0.4643, r, = 0.5403, r,a, 一 0.0760,r,，。 = 0.0269, ti zn E 
检验 : rs G= 1,3,4,6,7),P > 0.10;(7 = 2,5),0.05 < P < 0.10 表明 ,各 气象 因子 对 
午 指 数 的 相关 关系 不 显著 (P > 0.05)， 即 没有 一 种 因子 是 单独 、 显 著 地 起 作用 。 而 回归 
分 析 结 果 (7) 表 明 ,，7 个 气象 因子 对 重 数 量 有 显著 影响 (CP < 0.01), 是 共同 起 作用 的 。 由 
拟 合 值 表 3， 重 数量 残 差 均 在 y,— $4€1—26, 26] = [—0.7628, 0.7628] 置信 和 带 内 
(P < —26 < y, — J} < 26) 一 0.95，6 为 标准 误差 ), 表 明 气 象 因子 能 够 较 好 地 拟 合 草 指 
数 。 
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7 个 气象 因子 对 了 蚤 指数 能 够 综合 起 作用 CP < 0.01) 并 不 等 于 同等 看 符 一 切 因 子 。 无 
论 是 从 最 优 回归 子 集 还 是 逐步 回归 , 入 选 气象 因子 均 为 相对 湿度 x;:、 地 表 温 度 x, 和 气压 
x. 由 于 气象 因子 单位 不 一 致 ,不 能 对 系数 进行 比较 , 考虑 做 标准 化 变换 后 所 求 的 标准 回 
归 方 程 (8), 结 果 为 相对 湿度 和 地 表 温 度 的 系数 分 别 为 0.7154 和 0.6388, 是 接近 于 气压 系 
数 0.3610 的 2 倍 ; 又 在 全 部 回归 得 选中， 相对 湿度 和 地 表 温度 人 选 次 数 依 次 居 第 一 、 二 
位 ,表明 对 了 蚤 指数 影响 最 大 的 气象 因子 是 相对 强度 和 地 表 刘 度 。 

从 表 2 看 ,降水 对 重 指 数 的 影响 较 少 ,但 降水 是 影响 长 爪 沙 鼠 的 重要 气候 条 件 ” SERT 
相对 强度 、 党 内 齐 度 都 受降 水 影响 , 故 也 可 影响 荔 数量。 

光 是 一 个 重要 的 生态 因子 , 它 有 多 方面 的 作用 ,其 中 以 太阳 光 为 最 重要 。 由 于 跳蚤 怕 
光 , 喜 欢 钻 在 宿主 的 毛 层 中 , 以 避免 光线 刺激 ,因此 ,反映 在 回归 模型 中 所 起 的 作用 较 小 。 

跳 重 数量 的 变动 , 受 多 种 因素 的 影 秽 。 除 各 方面 都 依赖 于 宿主 外 ,气象 因子 x 一 x; 对 
KERALAN (P< 0.01)。 但 仅 由 气象 站 所 提供 的 数据 是 不 足以 解决 重 数 量 的 
波动 下。 利用 线性 模型 ,从 数量 上 研究 了 对 重 指数 影 啊 的 一 个 方面 ,而 现实 中 的 某 些 气 
象 因子 对 其 影 啊 可 能 是 非 线性 的 。 如 1984 年 的 拟 合 值 残 差 落 在 置信 带 内 是 满足 要 求 的 ， 
但 拟 合 值 是 负 值 就 是 一 个 不 合理 现象 , 它 表 明 可 能 有 某 种 非 线性 关系 。 
2.3 在 将 宿主 和 寄生 因子 、 气 象 因子 作为 两 组 变量 研究 时 ，(〈1) 由 于 典型 相关 变量 w 与 
的 关系 密切 (P < 0.001), 因 中 鼠 重 系数 的 比例 为 0.3174:0.7608 王 1:2.4， 知 气象 因 
子 对 重 指 数 的 密切 天 系 优 于 鼠 密 度 。 在 rn 中 ,起 作用 较 大 的 依次 为 蒸发 量 、 相 对 湿度、 地 
表 温 度 、 气 温 等 , 即 温 湿 因子 显著 影响 重 指 数 。 菊 发 量 在 vi 中 起 的 作用 大 , 其 原因 可 能 是 
在 鄂尔多斯 范 漠 草原 ,由 于 气 修 干燥 引起 ,但 集 内 温度 不 会 明显 增高 ， 这 样 可 能 对 么 繁殖 
有 力 。(2) 典 型 相关 变量 w 与 v; 的 相关 检验 已 < 0.25, 从 vz; 看， 起 作用 的 主要 为 降水 ， 
又 降水 与 w; 中 鼠 密 度 同 号 ,表明 降水 增多 时 ,局 密度 上 升 ， 反 之 亦 然 。 注意 到 Ux 指 
数 为 负 号 , 即 当 鼠 数量 增加 时 , 重 指数 可 能 下 降 。 这 是 因为 ,在 长 爪 沙 鼠 台 增 年 份 ,由 于 降 
KRZ, 洞 内 湿度 太 高 , 重 死 亡 数量 太 多 ,繁殖 数量 大 大 降低 , 鼠 与 入 数 量 可 能 呈 负 相关 趋 
Bo 
24 重 指 数 与 宿主 数量 、 气 象 因子 的 关系 密切 ， 它 还 与 动物 鼠疫 流行 与 否 关 系 密 切 
(P < 0.01)。 这 是 因为 ，1975、1987 一 1989 年 为 长 爪 沙 鼠 动物 鼠疫 流行 年 份 。 将 动物 鼠 
疫 流 行 年 和 未 疲 行 年 数量 化 为 1 和 0 , 计算 重 指 数 与 相应 年 份 间 的 相关 系数 > 一 0.6876 
> ros (11) 一 0.6835。 因 此 ,进行 卫 数 量 动态 研究 ,对 流行 病 学 分 析 和 鼠疫 动物 病 预 测 预 
报 ,同样 具有 重要 意义 。 
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Abstract According to the population of Meriones wunguiculatus, flea index and 
data of seven meteorological factors in Etuoke Banner and Etuokeqian Banner, Inner 
Mongolia Autonomous Region in 1975—1989, models of linear and curve regression 
between the flea index and the population of M. unguiculatus were obtained respec. 
tively. The optimum regression subsets of multipic linear regression were conducted 
between the flea index and meteorological factors. And then, the model of standard 
regression was also obtained. The canonical correlation analysis explained such rela. 
tions of the flea index and their hosts population to the meteorological factors. The 
conclusions are as follows: The population of M. wnguiculatus was the main factor 
affecting flea index. The meteorological factors exerted summing.-up influence on the 
flea population, and yearly average relative humidity and temperature in the field 
were the two main factors with influence. Moreover, the flea index was influenced 
by the meteorological data greater than the population density of their hosts. 
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